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We study the properties of a spin-density-wave antiferromagnetic mean-field ground state with 
d-wave superconducting (DSC) correlations. This ground state always gains energy by Cooper 
pairing. It would fail to superconduct at half-filling due to the antiferromagnetic gap although its 
particle-like excitations would be Bogolyubov-BCS quasiparticles consisting of coherent mixtures of 
electrons and holes. More interesting and relevant to the superconducting cuprates is the case when 
antiferromagnetic order is turned on weakly on top of the superconductivity. This would correspond 
to the onset of antiferromagnetism at a critical doping. In such a case a small gap proportional to 
the weak antiferromagnetic gap opens up for nodal quasiparticles, and the quasiparticle peak would 
be discernible. We evaluate numerically the absorption by nodal quasiparticles and the local density 
of states for several ground states with antiferromagnetic and d-wave superconducting correlations. 
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I. INTRODUCTION 

Ever since the discovery of high temperature 
superconductivitjii it was proposed that the supercon- 
ducting correlations might already exists in the antifer- 
romagnetic Mott insulator^. The origin of the supercon- 
ducting correlations was ascribed to the large Coulombic 
interactions in the undopcd materials. The only other 
large energy scale in the materials is phononic"^. 

While the microscopic origin of superconductivity re- 
mains a matter of debat8^*^*Sil, there is growing experi- 
mental evidence that the quasiparticles are Bogolyubov- 
BCS quasiparticles. Bending back of photoemis- 
sion bands^ fits quantitatively the BCS-Bogolyubov 
modelSiii. Scanning tunneling microscopy finds coher- 
ent quasiparticles that disperse as a coherent mixture 
of particles and holeaii*iS. The particle and hole am- 
plitudes in these experiments and in inverse photoemis- 
sion experimentsiSiH fit accurately to the theoretical 
Bogolyubov-BCS values calculated from the dispersion 
and gap measured in the normal and superconducting 
materials respectively. 

Regardless of whether the origin of superconducting 
correlations is exotic Coulombic physics or some more 
conventional mechanism, it is clear that the cuprates are 
BCS paired superconductors. This does not mean that 
the Coulomb interactions do not matter. Rather, the in- 
teresting and contradictory physics for underdoped ma- 
terials is the result of Coulomb degradation of the super- 
fluid densitySiii^ and order parameter competition be- 
tween superconductivity and correlated electron ground 
statesi^ii^. The degradation of the superfluid density 
leads to suppressed Tc due to a phase instability of the 
superconducting order parameteri2iiSiiLi2iiS. There are 
several Coulomb stabilized competing ground states such 
as orbital antiferromagnetism^^, stripe or charge density 
wave ground states^*^ and perhaps electronic liquid crys- 
tal phases^-. Regardless of which of these competing 
ground states are realized, there is strong experimental 
evidence for incommensurate electronic orderingSSi, either 
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FIG. 1: Spectral Function for the SDW ground state without 
(la) and with (lb) superconducting correlations built in. 



static or incipient. The evidence seems more consistent 
with charge-density-wave or stripe order. 

In the present work we will study the physics of an 
antiferromagnet with a strong d-wave Cooper pairing in- 
teraction. We do not speculate as to the origin of this 
superconducting interaction except to point out that in 



2 




FIG. 2: In the upper left side we draw two bands separated 
by a gap, with the lower band partially filled appropriate to 
a metal or superconductor. In the upper right side we draw 
the situation encountered for the electron fluid under the ac- 
tion of an electric field. The lower sketch illustrates the sit- 
uation appropriate to an insulator, where the lower band is 
completely filled, making conduction impossible regardless of 
Cooper correlations. 

such a model it competes with the Coulombic antiferro- 
magnetic physics. Both the superconductor and the anti- 
ferromagnet are studied in the mean field approximation. 
While one can doubt the validity of such an approxima- 
tion at a phase transition point, it will be qualitatively 
correct within the ordered phases. 

Cooper pairing leading to a BCS ground state is an 
instability of a Fermi liquid ground state. In this study 
we apply the BCS approximation to a spin density wave 
(SDW) insulating ground state as it exists in the cuprates 
at half filling. The resulting ground state has Cooper 
pairing yet it fails to superconduct due to the SDW insu- 
lating gap. Next we will review some well known facts in 
order to understand how a state with Cooper pairs does 
not superconduct. Before doing so we emphasize that 
this only happens as a consequence of having a completely 
filled insulating band. 

When an electric field is applied to a metal, it con- 
ducts dissipatively. The way this happens is that the 
center of mass of the Fermi sea gets displaced upward 
in the unfilled metallic bandS^. Ohmic dissipation oc- 
curs because newly filled electronic states at the top of 
the Fermi sea get scattered into newly empty electronic 
states at the bottom of the Fermi sea due to the lack of 
rigidity of the Fermi liquid ground state (see figure 2). 
When there are Cooper correlations, the electron liquid 
gets displaced upward in the band too, but as long as 
the displacement in energy within the band is less than 
the superconducting gap. Cooper pair correlations make 
the electron liquid rigid, thus preventing scattering and 
dissipation. For the case of an SDW ground state at half 
filling with Cooper pairing correlations there is no super- 



conductivity as the electron fiuid cannot move upward in 
the band for the band is full and there are no electronic 
states to be filled unless one excites across the insulating 
gap and into the conduction band (see figure 2). 

That the SDW insulating ground state with d-wave 
pairing interactions has Cooper pairing in the ground 
state can be seen from figure 1. In figure la we plot 
the spectral function for the SDW ground state with no 
superconducting correlations. In figure lb we plot the 
spectral function for the SDW ground state with super- 
conducting correlations. In the ground state with both 
superconductivity and antiferromagnetism, the separa- 
tion between the coherence peaks is bigger as it gets con- 
tributions from both the SDW and superconducting gap. 
A prediction of this model is that the quasiparticles will 
be coherent with an electron and a hole component in 
agreement with the BCS-Bogolyubov model. 

The SDW ground state with d-wave Cooper pairing 
(SDW-DSC) will become superconducting when doped. 
At the mean-field level, without worrying about self con- 
sistency, the chemical potential will jump to the ap- 
propriate band and there will be a low superfiuid den- 
sity superconductor. Whether this physics is correct 
for the cuprates is controversial. There is experimen- 
tal evidence for the chemical potential staying pinned at 
midgap due to spectral redistribution of states toward 
midgap states24. There is also experimental evidence for 
chemical potential shifts in the cuprates, in the same 
way as in regular semiconductor materials2^. Indepen- 
dently of whether the SDW-DSC ground state has chem- 
ical potential shifts or not, the physics of an insulator 
with Cooper pairing correlations is interesting. For our 
study we have the cuprates in mind. For these materials, 
some phenomenology of this form seems to apply^, but 
it would be interesting if this physics were to be realized 
in nature irrespective of the cuprate problem. 

In the present work we will flip the problem around. 
We will start with a d-wave superconductor (DSC) and 
begin turning on SDW antiferromagnetic order on top of 
the superconductivity. In this limit, the complications 
mentioned in the previous paragraph are nonexistent. A 
slow turning on of SDW order on top of the supercon- 
ductivity will show up as a shift of the antinodal gap and 
a gapping of the nodal quasiparticles. The latter should 
be a signal much easier to pick out than the gap shift. 
The gapping of the nodal quasiparticles is not a unique 
prediction of antiferromagnetic ordering on top of the su- 
perconductivity, as such a gapping can be produced by 
disorder. On the other hand, the coherence of the gapped 
"nodal" quasiparticles would be nonexistent for a disor- 
dered gap and is thus a unique signature of antiferromag- 
netic ordering developing on top of the superconductivity. 
Therefore, if a quasiparticle peak is discernible, and the 
broadening is less than the disorder-induced broadening 
(^ jlmt^x^ ~ 350 meV for Ax ~ 1 nm, appropriate 
to the cuprates), then the gap is a long range ordered 
gap and not a disordered gap. Another unique signature 
of an SDW gap is that the gap will open exactly at the 
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FIG. 3: Gapping of the nodal quasiparticles pole as the SDW 
order develops. 



doping where the antiferromagnetism starts. 
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There are experimental suggestions of antiferromag- 
netism competing with superconductivity in the deep un- 
derdoped regime in the cuprates. For example, measure- 
ments show the nodal quasiparticle peaks surviving right 
up to the doping where antiferromagnetism starts. The 
spectral weight of such peaks diminishes with decreas- 



ing doping, consistent with spectral weight being robbed 
from the superconducting long range order by a compet- 
ing long range order such as antiferromagnetism^. If one 
looks in the antiferromagnetically ordered dopings, there 
are experimental suggestions of a competing order pa- 
rameter that conducts efficiently. Most strikingly, there 
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FIG. 4: Shift of the antinodal gap as the SDW order develops. 



are nieasurements of metallic conduction even below the 
Neel ordering temperature^i. 

The gapping of the nodal quasiparticles pole as the 
SDW order develops on top of the superconductivity is 



shown in figure 3 for different values of the SDW gap. 
The reason we only have a quasiparticle sharp pole is that 
we have not modeled the realistic electronic self energies 
relevant to the cuprates as they are irrelevant to the point 



of principle we are making. Their only effect will be to 
broaden the quasiparticle peaks and add an incoherent 
background with the phenomenological features. In fig- 
ure 4 we plot the shift of the antinodal gap as the SDW 
gap turns on. In figure 5 we plot the spectral density 
of states in a d-wave superconductor as the SDW gap is 
turned on. The superconductor with no SDW gap does 
not have a true gap because of its d-wave symmetry. This 
is seen in the familiar V-shaped collapse at zero energy. 
As the SDW gap is turned on, we see the V-shape flatten 
and expand as a signature of the opening of the antifer- 
romagnetic gap. 



II. HUBBARD MODEL WITH D-WAVE 
ATTRACTIVE INTERACTIONS 

For the cuprate problem, the two large effects are 
the antiferromagnetic, or Coulombic, physics and the 
strong superconductivity. Hence we will start from a phe- 
nomenological Hamiltonian which is a Hubbard model 
with a d-wave electronic interaction. This interaction 
will give rise to d-wave superconductivity when we make 
the mean-field BCS approximation. The Hamiltonian is 



/-^^ ' ' fe,(T fe,"" _/V ^ '=liT '^l+^'T k2+q,i *:2,i 

k.(7 ki.k2,q 

ki ,k2 
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FIG. 5: Spectral density of states for a d-wave superconductor 
as the SDW gap increases. 
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where , cr are the electronic creation and destruc- 

tion operators with momentum k and spin a, eg is the 
kinetic energy, fi the chemical potential, and U is the 
Hubbard repulsion. We are working in a spatial lattice 
with TV sites. The last term is an electronic interaction 
chosen in the reduced BCS formS, which will be used 
to stabilize superconductivity. In order to have d-wave 



superconductivity we choose V{ki,k2) = Vq {coskix — 
coskiy) (cosk2x — cosk2y) ■ This phenomenological Hamil- 
tonian can have a mean-field SDW ground state and a 
mean-field DSC ground state. It can be used to study 
the turning on of DSC correlations on top of an SDW 
ground state, or the turning on of SDW order on top of 
the superconductivity. 
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We will analyze this Hamiltonian by imposing an SDW 
mean field condition, which is stabilized by the Hubbard 
term. This will be followed by a DSC mean-field con- 
dition, which is stabilized by the reduced BCS d-wave 
interaction. While the use of two mean-field conditions 
is not common, it has important precedents. It was used 
by P. W. Anderson-— in his study of the role of plasmons 
in restoring gauge invariancc to the BCS ground state. In 
this work he invented the Anderson-Higgs mechanismSS. 
He solved for the properties of the electron system im- 
posing a mean-field condition on the electron density, as 
in the study of electron correlations by Sawada, et aS^ 
and a BCS electron pairing mean-field conditionsli. 

The Hubbard interaction stabilizes the mean-field or- 
der 



/-^^ k+O 



Cr 

k+Q,a k,, 



(2) 



where Q = (tt, tt) is the commensurate ordering wave 
vector and S is the average magnetic moment per site. 
Other ordering wave vectors are possible for spin and/or 
charge, i.e. stripe, order parameters but we do not con- 
sider them in our study. When we impose this condition 
on the Hamiltonian and neglect fluctuation terms, the 
Hamiltonian becomes 



III. BOGOLYUBOV DIAGONALIZATION OF 
THE MEAN-FIELD HAMILTONIAN 

We know diagonalize the Hamiltonian by the Bo- 
golyubov methodiSi. We will do this in two steps. First 
we diagonalize the SDW part. In order to do this more 
conveniently, we will split the momentum sums into sums 
over the reduced magnetic zone. The Hamiltonian is then 



+ UNS^-^+y'Aj: (ci ct- 



^k+Q,a'^k+Q,a 



k,<^ k+Q,a k+Q,ij' 

-2aUSci - cr ] 



^-k.l'^k.l 



(6) 



where the prime on the summation sign means that the 
sum is restricted to the wave vectors in the magnetic 
zone, = (eg + e^^^)/2 and eZ = (eg - e^^Q)/2. The 
last term in the superconducting interaction is negative 
because Ag^^ = — Ag. In order to diagonahze the mag- 
netic part we define the Bogolyubov operators 



Ti: = y (er - ^) ct Cr: + UNS^ - U Secret - c,- 



k,a 



k.a 



ki ,k2 



^ fe2,i%:T 



(3) 



We see that by ordering antiferromagnetically we gain 
variational energy —UNS'^ if self-consistency can be 
achieved. We next impose the mean-field d-wave Cooper 
pairing 



bt — oitcr 

k,<7 k k,< 



""PkH+Q^a 



^k+Q,rj ~ ^k^k+Q.a + '^Pk^k,a 



If we choose 



(7) 



(8) 



-I 



Ag^ = {cosk2x - cosA:2y)Vb^(cosfci:c - coskiy){c^^^ ^c^^^ ^ 

(4) 



fci 



= /^n{cosk2x ~ COsk2y) 
Then the Hamiltonian becomes 



k,tT^k,, 



UNS^ 



USy^aci - c.- 

Z-^ k+Q,a k,i 



_^ + y-A. (ct cK 



^-k,l'^k,T 



(5) 
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the Hamiltonian becomes 



k.a 



u)(bi bt +hi - 6.- ^ ) 

f^l^ i. „ k.a k+Q,a k+Q,iyl 



(10) 



kM k,a 

Er(hi bt -hi - 6,-^^ )\ + UNS^ 

k,a k,(T k+Q,a k+Q,a) \ 
^ + y'Ar (bl bK +b t,bt, 



We see that if the phenomenological d-wave interaction 
is attractive, i.e. Vg < 0, we gain variational energy 
Aq/Vq by Cooper pairing regardless of whether we have 
ordered antiferromagnetically or not. Of course, if we are 
at half filling, the material will be insulating irrespective 
of the presence of Cooper pairs, as we would have to 
excite quasiparticles across the SDW insulating gap in 
order for conduction to take place. 



5t . &t_ 
fe+Q,T -k-Q,l 



b_k.Q,ih+Q,i) (11) 



Our last step to diagonalize the full Hamiltonian is the 
Bogolyubov diagonalization of the leftover superconduct- 
ing part by defining the canonical operators 



Br 



U+br 
k k,<j 



av+bi 

k k,a 



(12) 
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If we choose 



k+Q,a 




k k+Q,g- 



(13) 



(14) 



how they rob spectral weight from each other. In the 
present work we do not worry about self-consistency but 
this is very important and interesting. Hence it will be 
part of future work. In this work we will only concentrate 
on the spectral properties of the system with coexisting 
SDW and DSC correlations. 



IV. GREEN'S FUNCTION FOR THE SDW-DSC 
HAMILTONIAN 




the Hamiltonian then becomes 



(15) In the present section we will write down the expres- 

sions for the Green's functions for a system with SDW 
order and d-wave Cooper pairing in each of the SDW 

, . bands. The expression for the retarded Green's function 

^ ' or the propagator is^i 



n 



y'\E+Bi 

L k k,< 



Br +EZBi - Br^n 

t','^ k k+Q,a k+Q.a 



A2 

-I- UNS^ ~-rr+ constants (17) 

We see that we have two separate superconducting bands 
with dispersions Et and EZ . This happens because the 

k k 

SDW ordering has split the nonintcracting band, i.e. the 
system with U — 0. Of course, if the SDW gap were to 
collapse, the two bands would merge into one supercon- 
ducting band. If we look at it from the opposite perspec- 
tive, we see that when we turn on the SDW order, there 
will be an insulating gap. We have shown in figure 2 how 
this gap opens up at the node as calculated in the next 
section. 

The Bogolyubov transformations can of course be in- 
verted to yield the electron creation and destruction op- 
erators in term of the Bogolyubov eigenoperators (Bj:) 
of the system. From this we can evaluate the self consis- 
tency or "gap" equations (jH) and Q. We obtain 



G{x,x , t) 



-iE„t/h 



^(-t)e'^"*/''(V'ol4^JVn)(^n|c.-.|^o) } (20) 



2m 



duj 



where 



1 



(21) 



k 



(22) 



n labels the eigenstates of the system with energies En 
and the ground state energy i^o has been chosen to be 0. 

From the time Fourier transform of the Green's func- 
tion above we obtain the local retarded propagator in the 
energy representation 




from the antiferromagnetic self-consistency condition 
and 



+ \2 



irj 



E-Et ^ 

k 



E + Et -ir] E + E^ 



IT] 



E-EZ 



irj 
(23) 



2 



k 



, ^2 / 1 1 

cosky) — X H 



(19) 



from the superconducting sel-consistency or gap equa- 
tion Q). The negative sign on the left is consistent with 
Vb < as is necessary to stabilize superconductivity. We 
see that these are two coupled equations for the super- 
conducting gap parameter Ag and the spin moment mag- 
nitude S in the antiferromagnet. Their solution will con- 
tain information about how the two orders compete and 



The local spectral density function follows fromai 



A{x, E) = ImG(f , X, E) 

TT 



(24) 



All of our density of states are calculated from these ex- 
pressions in a 1000 x 1000 momentum lattice with an 
energy resolution of 0.01. We choose the hopping energy 
scale to be 1, so all energies are measured in hopping 
units. When we have superconductivity we choose the 
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gap to be 0.3. The antiferromagnetic gap is chosen any- 
where between and 0.6, usually with jumps of 0.1. We 
have nearest neighbor hopping only. These values need 
not be realistic; they are just chosen to illustrate the ef- 
fect. 

Similarly, if we Fourier transform the Green's function 
(|20|l in both time and space, we obtain the retarded prop- 
agator in the wavevector energy representation. 



G{k,E) = 



IT] 



E-EZ 

k 



17] 



E + Et -ii-j E + E^ 



(25) 



From this formula we calculate the absorption strength 
vs. energy for the nodal quasiparticles. We do this by 
simply fixing k to be at the node and plotting the spectral 
densitv^ 



A{k, E) = ImG'(fc, E) 

TT 



(26) 



vs. energy. Energy units, values and uncertainties are 
chosen as described for the local density of states. 



V. CONCLUSIONS 

We studied a mean field Hamiltonian with two mean 
field order parameters. The Hamiltonian contains a spin- 
density- wave antiferromagnetic mean field stabilized by a 
Hubbard interaction and a d-wave Cooper pairing mean 
field stabilized by a phenomenological d-wave interac- 
tion. The two order parameters can coexist and the SDW 
ground state always gains energy by Cooper pairing when 
the d-wave interaction is attractive and nonzero. The 
SDW ground state with Cooper pairing fails to super- 
conduct at half-filling due to the antiferromagnetic gap. 
Its particle-like excitations are Bogolyubov-BCS quasi- 
particles consisting of coherent mixtures of electrons and 
holes. 



Of greater interest and relevance to the superconduct- 
ing cuprates is the case when antiferromagnetic order is 
turned on weakly on top of the superconductivity. This 
would correspond to the onset of antiferromagnetism at 
a critical doping. In such a case a small gap proportional 
to the weak antiferromagnetic gap opens up for nodal 
quasiparticles, and the quasiparticle peak would be dis- 
cernible. While the gapping of the nodal quasiparticle 
could be caused by a large enough disorder, such a dis- 
order would broaden the quasiparticle peak so much as 
to make it invisible. A unique signature of antiferromag- 
netic gapping of the nodal quasiparticles is that it will 
turn on always at the doping when antiferromagnetism 
starts while disorder gapping will turn on at different 
sample dependent dopings. 

We wrote down the exact expressions for the Green's 
fimction for the system with coexisting SDW and DSC 
order parameters. These are evaluated numerically in a 
1000 X 1000 momentum lattice with .01 energy resolu- 
tion in units of the lattice hopping. From the imaginary 
parts of the Green's functions we obtained the absorption 
by nodal quasiparticles and the local density of states. 



In our work we did not worry about having self- 
consistency. This neglect does not affect our results when 
the two order parameters are nonzero, but it will affect 
whether the order parameters are nonzero or not, and 
what the gap values are. Self-consistency will be im- 
portant in studying how the two order parameters com- 
pete and if and how they steal spectral weight from each 
other. Self-consistency might also affect how the SDW- 
DSC ground state behaves when doped from half-filling. 
Intuitively one expects chemical potential shifts, but it 
is not certain that this would be the case. All these is- 
sues should be studied carefully and we postpone them 
for future work. 
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